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We consider defect mediated melting of Fulde-Ferrell-Larkin-Ovchinnikov (FFLO) and pair den- 
sity wave (PDW) phases in two dimensions. Examining mean-field ground states in which the 
spatial oscillations of the FFLO/PDW superfluid order parameter exhibit hexagonal lattice sym- 
metry, we find that thermal melting leads to a variety of novel phases. We find that a spatially 
f"^ , homogeneous charge six superfluid can arise from melting a hexagonal vortex-anitvortex lattice 

£\j ■ FFLO/PDW phase. The charge six superfluid has an order parameter corresponding to a bound 

state of six fermions. We further find that a hexagonal vortex-free FFLO/PDW phase can melt 
to yield a conventional (charge two) homogeneous superfluid. A key role is played by topological 
defects that combine fractional vortices of the superfluid order and fractional dislocations of the 
lattice order. 
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The interplay between solid order and superconducting/superfiuid order has become an issue of tremendous in- 
terest in a variety of physical systems. The putative supersolid phase of 4 Hei has provided a strong motivation to 
understand the relationship between these two orders. Additionally, the Fulde-Ferrell-Larkin-Ovchinnikov (FFLO) 
■ phased, recently observed in ultracold 6 Li atom systems^, provides another compelling example. This phase exhibits 
translational symmetry breaking through the formation of a paired fermion superfluid lattice over which the spatial 
average of the superfluid order is zero. Related pair density wave (PDW) states, a generalization of FFLO phases, 
are relevant in CeCoIng£~— and in the cuprates Lai_ a: Ba a ;Cu04 and Lai.g-xNdo^Sr^CuO^^. Also related to both 
PDW and FFLO phases are vortex-antivortex (v-av) lattice phases. Such phases have been discussed in the context 
in 2D superfluid 4 Heii and superconducting thin films 12 , where, at high temperatures there exist thermally excited 
vortices and antivortices. It is argued that if the density of these vortices is sufficiently high then a v-av solid phase 
will appear rather than a v-av liquid. Recently, a staggered vortex phase (a specific type of v-av phase) has been 
observed in an optical square lattic e) 13 ' 14 . This diverse set of physical systems underlies the ubiquity and importance 
£j • of the interplay between these orders. 

Very recently, a quasi-one dimensional FFLO phase has been engineered in cold atomic gases^. This system allows 
for the opportunity to examine recent theoretical predictions of the Larkin Ovchinnikov (LO) phase, in which the 
mean field superconducting order parameter breaks translational symmetry in one direction (for example, the pairing 
gap A(x) — A cos(qx)). In particular, recent theoretical work has focussed on the consequences of the U(l) x U(l) 
symmetry that exists in the free energy due to translational and gauge invar ianc o 1 ' 1 5 ~ — . This symmetry implies 
\ the existence of fractional vortices in addition to the usual vortices and dislocations that would be anticipated. The 
fractional vortices have a superfluid phase winding of 7T, 1/2 the usual vortex phase windin g 10 1 15 ~ — . These fractional 
vortices are thus called 1/2 vortices. This 1/2 phase winding is accompanied by a 1/2 dislocation, which leads to 
an additional sign change in the order parameter, so that the order parameter remains single valued when the 1/2 
vortex is encircled. In two-dimensions (2D), the existence of these 1/2 vortices are closely related to the existence 
of a spatially homogeneous superfluid condensate that corresponds to a bound state of four fermions: a charge 4e 
superflui d 10 ' 15 (in rotationally invariant superfiuids, this phase always appears in two and three dimensions 1 ^). 
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The recent discovery of a LO phase in quasi-lD cold atoms systems provides an ideal opportunity to examine the 
physics discussed above. It further indicates that that FFLO phases confined to 2D are also likely to be realized in 
cold atoms. Mean field theories of the FFLO phase in 2D predict not only an LO phase, but a variety of other stable 
FFLO phases^. In many of these phases, superfluidity is spatially modulated with an underlying hexagonal or square 
lattice. These results lead to some more general questions about PDW phases: Are there other possible fractional 
vortices? Are there other exotic phases? Here we examine two such 2D FFLO phases with underlying hexagonal 
lattices that have been found in microscopic theories 1 ^. The theory we develop is relevant not only to FFLO phases, 
but also to PDW and v-av lattice phases discussed above. Our most interesting results are: the existence of a spatially 
uniform charge six superfluid phase, in which quasi-long range order appears only in an order parameter corresponding 
to a bound state of six fermions (this results from a mean-field v-av lattice phase in which 1/3 vortices exist); and 
the existence of spatially uniform charge two superfluid phases (this results from a mean-field phase in which there 
are no fractional vortices). These phases, in addition to the charge four superfluid 10 - 5 and non-superfluid density 
wav e 10 ' 15 ' 17 found earlier in stripe-like PDW and FFLO phases indicate that the physics of such phases is much richer 
than previously anticipated and offer the possibility to see never before seen states of matter. 
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FIG. 1. Directions of Q t used in the text. 



I. GINZBURG LANDAU WILSON THEORY 



We consider two symmetry groups for the normal state which will serve to define the FFLO/PDW order. Both are 
two dimensional (2D): the first is an isotropic normal state with cylindrical rotational symmetry and the second is the 
2D space group P6m, the group of a triangular lattice. For simplicity, our development and emphasis will be on the 
group P6m and we will state results for the case with cylindrical symmetry. For FFLO/PDW order appearing at a 
wavevector Q, the order parameter is defined by the irreducible representations of Gq (the set of rotation elements that 
conserve Q) and the star of the wavevector Q in P6m (set of wavevectors symmetrically equivalent to Q)2£. We choose 
the wavevector Q 2 = if"^f (0> 1) (where a is FFLO/PDW lattice constant), which is invariant under the rotation group 
Gq = {E, C2 y , cr z ,a x } with Ciy the 180°-rotation around the axis (1, 0), u 2 and a x the mirror operations perpendicular 
to the 2D plane and the plane perpendicular to (1,0), respectively. The irreducible representations of Gq are all one- 
dimensional. The only situation in which the particular irreducible representation of Gq is relevant is when there is 
a spatially uniform Q = supcrfluid order also present (for example, corresponding to usual Cooper pairs)^. This 
situation can be accounted for easily and, for this reason, we consider explicitly the identity representation in the 
following (for which the order parameter is unchanged under the action of any clement of Gq). To define the additional 
order parameter components at the wavevectors in the star of Q we use the elements {E, Ce, Cf , Cf , C|, Cf }, these 
give the star of Q 2 , {£?2! — Q\i Q3; — Q2 Qii as shown in Fig. [1] This then defines a superconducting order 

parameter with six complex components which we define as A = (Aq x , Aq 2 , Aq 3 , A_q x , A_q 2 , A_q 3 ). We take 
Q x = ^L-2=(y / 3/2, -1/2), Q 2 = 2^-^(0,1), and Q 3 = Q, - Q 2 so that the superfluid order is unchanged by 

the translations cti = a(l,0) and a 2 = a(l/2, y/3/2) (note that these are not translation vectors of the underlying 
microscopic triangular lattice). We consider the case that these translations vectors are not commensurate with the 
translation vectors of the underlying microscopic triangular lattice (as is the usual case with FFLO phases). With 
these definitions, the symmetry properties of the order parameter are given as follows: under a microscopic translation 

T, A Qj -> e^J'^AQj ( Ag^ -» e^^j'^Ag.) and under a time-reversal operation Aq^ -> A*_q.- Moreover, under 
point group symmetries we have that (Aqj , Aq 2 , Aq 3 , A_Qj , A_q 2 , A_q 3 ) transforms to 



The GLW free energy is constructed by requiring invariance under the group P6m, U(l) gauge symmetry (under 
which Aq ; — > e Ag i ), and time- reversal symmetry. These symmetry operations are given above and the resulting 
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(A-Q 3 , A_ Ql , A_q 2 , Aq 3 , Aq 1 , Aq 2 ) 

( A Qi 1 A Q 2 5 Aq 3 , A_ Ql , A_q 2 , A_q 3 ) 

( A Q 3 > A Q 2 j A Qi 7 A -Q 3 ) A -C?2 > A -Qi ) 
(Aq 3 , Aq 2 , Aq 1 , A_ Q3 , A_ Q2 , A_ Ql ) 
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Phase 


(Aqj , Aq 2 , Aq 3 , A_ Ql , A_q 2 , A_q 3 ) 


Free Energy j3 


*LO 
*2Q 
^ v — av 
*3Q 
*4Q 
*A 

#6Q,1 

*hc 
*hc,2 
*6Q,2 


e i8 (l, 0,0, 0,0,0) 
^(e^,0,0,e-^\0,0) 
^(e^,0,0,0, e -^,0) 
^( e **i )e *<h }e -H<h.+M } 0)0 ,0) 

^ (e^ 1 , ie"* 2 , 0, e"^ 1 , ie"^ 2 , 0) 
e' 8 ( e i » >1 ,e i<# ' 2 ,ae-* Wl+ ' >2) ,fee- i,# ' 1 ,6e- i '* 2 ,ce lWl+02) ) 

y 2+a 2 +2b 2 +c 2 

^(ie'^ 1 e i7T ^ 3 e ir ^ 2 e~ i7V /3 e -i(.<t>-t+4>2) _^ e -i<t>i e W3 e -s<fe e -W 3 giW>i+'fe)) 
e"> (ie i01 ,e^ 2 ,ae- 4Wl+02) , , fee""* 2 , ce iWl+te) ) 

V2+2b 2 + a 2+ c 2 v ; 


Pi 

ft + ft/4 
Pi +/? 4 /4 
ft + ft/3 
Pi - 4 W f- Wa \ (P3<0,2\p 4 \ <ft) 

ft + (ft +ft +ft-ft)/8 

pi + fa/12 + (ft + Pi + ft)/6 (7<0) 
ft + (ft-ft)/12 + (ft + ft)/6 ( 7 <0) 
no analytic solution found (ft < 0, 7 < 0) 

ft + ft/12 + (ft + ft + ft)/6 (7 > 0) 
ft + (ft - ft)/12 + (ft + ft)/6 ( 7 >0) 
no analytic solution found (ft < 0, 7 > 0) 







TABLE I. Possible FFLO/PDW ground states and associated free energy. The free energy is given by — a 2 /4ft The conditions 
in the brackets are necessary (but not sufficient) for the phase to exist. The phase factors 8, <f>\, and <f>2 are not determined by 
the free energy and lead to Goldstone modes of the FFLO/PDW phases. The parameters e, a, b, and c are determined by the 
free energy and are temperature dependent. 



free energy density is 

/ = -« E* I A Qi I 2 + ft (Ei I A Qi I 2 ) 2 + ft E, |Aq» I 2 1 a_ Qi i 2 

+ft(|A Ql | 2 |A Q2 | 2 + |A Ql | 2 |A Q3 | 2 + |A Q2 | 2 |Aq 3 | 2 + |A_ Ql | 2 |A_ Qa | 2 + |A_ Ql | 2 |A_ Q;! | 2 + |A_q 2 | 2 |A_q 3 | 2 ) 
+ft(|A Ql | 2 |A_ Q2 | 2 + |A Ql | 2 |A_ Q3 | 2 + |A Q2 | 2 |A_ Q3 | 2 + |A_ Ql | 2 |A Q2 1 2 + |A_ Ql | 2 |A Q3 1 2 + |A_ Qa | 2 |A Q3 1 2 ) 
+ft[A Ql A_ Ql (A Q2 A_ Q2 )* + A Ql A_ Ql (A Q3 A_ Q3 )* + A Q2 A_q 2 (Aq 3 A_ Q3 )* + c.c] + k x E, |VA Qi | 2 
+ K2 [^ 2 (|V + A Ql | 2 + |V+A_ Ql | 2 ) + (|V + A Q2 | 2 + |V + A_ Q2 | 2 ) + K|V+A Q3 | 2 + |V+A_q 3 | 2 ) + c.c.] 

(2) 

where c.c. means complex conjugate, V± = V x ± iV y , and v — e l27r / 3 . If the ground state solution has all six 
components unequal to zero, then Eq. [2] is not sufficient to completely specify the order parameter (there remains an 
unphysical U(l) symmetry in the solution). In this case the following free energy contribution is also required 

7 [Aq 1 A 02 Aq 3 (A_q 1 A_q 2 A_q 3 )* + A_ q1 A_q 2 A_ Q3 (A Qi A Q2 A Q3 )*]. (3) 

Unlike the case for tetragonal symmetry^ii, we are not able to analytically find the ground states of Eq.[2j However, 
through a combination of analytical and numerical analysis, we find that the Eq. [2] allows at least twelve possible 
global minima depending on the parameters ft . These are listed in Table U (note that these ground states also exist 
for a material that is cylindrically invariant). Of particular relevance are the states ^lo, and ^ v - av since these 
have all been found as stable ground states of cylindrically symmetric microscopic weak-coupling theories of the FFLO 
phased. The state ^lo has been previously studied and was discussed in the introduction. Consequently, in the 
following, we concentrate on the states \Pa and f„_„„. The state described by "J a is a superfluid triangular lattice 
(this state is closely related to ^hc, which is a superfluid honeycomb lattice without any vortices). The state *ff v - av 
is a v-av triangular lattice. These two states are depicted in Fig. [2] Prior to examining these two ground states in 
detail, we note that some of the other phases are also of interest. Perhaps the most interesting is ^kag- The local 
maxima of the superfluid density of ^kag form a Kagome lattice. Within the the hexagons of this Kagome lattice 
there are double vortices and there are single anti-vortices within the triangles of this Kagome lattice. The analysis 
that follows can be applied to any of the ground states listed in Table Q] 

II. SECONDARY ORDER PARAMETERS FOR THE PHASES AND * A 

In addition to the FFLO/PDW order parameters, there are secondary order parameters that play an important 
role in thermal melting and in distinguishing the different FFLO/PDW phases. In the mean field theory, these 
order parameters appear at the mean field phase transition in addition to the FFLO/PDW order. These secondary 
order parameters include density wave order, orbital angular momentum and spin density wave order (characterizing 
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A v-av 

FIG. 2. Hexagonal FFLO/PDW states considered in detail in this paper. The dots depict maxima in the magnitude of the 
superfluid density and the + (— ) symbols depict vortices of positive (negative) phase winding. Both states are stable 2D FFLO 
states in the weak-coupling theories with cylindrical symmetry. 



the vortex-antivortex lattice), and spatially uniform superfluid order. When thermal melting is considered, these 
secondary order parameters may become the primary order parameter and therefore play an important role in the 
theory. In the following, we characterize these secondary order parameters for the two states and ^ v - av in turn. 

The state is characterized by a spatially oscillating superfluid density with an underlying triangular lattice and 
a co-existing spatially uniform s-wave charge two superfluid order. Specifically, the secondary orders are: a spatially 
uniform conventional superfluid order tp s oc Aq i Aq 2 (A_q 3 )* + Aqj Aq 3 (A_q 2 )* + Aq 2 Aq 3 (A_q 1 )* and the density 
wave order pq 1 cx A_q 3 (Aq 2 )*, pq 2 oc A_q 1 (Aq 3 )*, and pq 3 oc A_q 2 (Aq 1 )*. The density wave order has the same 
lattice as the FFLO/PDW order. The appearance of a spatially uniform charge two superfluid is somewhat surprising 
for a FFLO/PDW phase. It is a consequence of the underlying hexagonal symmetry, it does not occur for the LO phase 
or for FFLO/PDW phases with an underlying square lattice. The existence of this order stems from the following 
coupling term in GLW free energy 

e{V,[A Ql (A_ Q2 A_ Qa )* + A Qa (A_ Ql A_Q 3 )* + Aq 3 (A_ Qi A_ Q2 )*] +c.c}, (4) 

since i/j s appears linearly, it must become non-zero at the mean- field FFLO/PDW transition. We note that, for the 
same reasons as 'Fa, the phases have ^hc,2 and ^kag have spatially uniform charge 2e p x + ip y and d x 2_ y 2 + id xy 
order, respectively, in addition to the FFLO/PDW order. 

The state describes a triangular v-av lattice. For ^ v - av the secondary order parameters are: a charge six 

superfluid order ipQ e oc Aq 1 Aq 2 Aq 3 and the orbital angular momentum (l z ) and density wave (p) orders il z ,K.i oc 
PKi oc Aq 2 (Aq 3 )*, il z ^ 2 oc p K2 oc Aq 3 (A Qi )*, andi£ z ,K 3 oc p K:5 oc Aq 2 (A Qi )* whereKi = Q 3 -Q 2 , K 2 = Qi~Q 3 , 
and K3 = Q 2 — Q 1 . The density wave order characterizes a hexagonal lattice that is rotated ir/2 and has a y3 shorter 
lattice vector than the FFLO/PDW lattice. The orbital angular momentum l z describes the v-av lattice that exists 
in this phase (the state ^a has no vortices, so that l z = 0). 



III. ELASTIC THEORY 



The undetermined phase factors 9, <fii, and 4> 2 in ^a and ^ v - av reveal an underlying £7(1) x £7(1) x £7(1) symmetry 
of the GLW free energy and correspond to elastic modes of these phases. Physically, this symmetry originates from 
the two phonon degrees of freedom u = (u x ,u y ) of the 2D FFLO/PDW lattice and the superfluid phase degree of 
freedom, 9. The phases 4>i and 4>2 are then <f>x = Qx ■ u and (j>2 = Q 2 ■ u. Under uniform phase shifts of 9 and u the 
order parameters transform as follows: 

A Q ^e i8+iQ - u A Q 
P Q ^e iQ u p Q 

I z ,k ^e lK u £ z , K . (5) 

At low temperatures, we can ignore fluctuations in the magnitude of the Aq^ and the effective Hamiltonian is governed 
by fluctuations in 9 and it. To lowest order, the relevant elastic Hamiltonians are given in Table [TTJ We have used 
the usual definitions ufj — ^(ff 1 + ^r) and u^- — ^(ff 1 — j^r)- Note that if cylindrical symmetry is assumed, 
then 7 = in the elastic Hamiltonians. However, if there is an underlying microscopic hexagonal lattice, then 7^0. 
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Phase 


Elastic Hamiltonian 


Vortex Charge 


Dislocation Charge 


^ v — av 


ip s (V0) 2 + f <4 + p(u%) 2 + 7 (u^) 2 
i Ps (V0) 2 + ±ul + M «) 2 + 7 («°) 2 


^ § d8 — |(m + n 2 + 7t 3 ) 


£ du — h(a, 0) 

I /„/« V3a\ 

= |(ni - ns ) 



TABLE II. Elastic Hamiltonians and the vortex and dislocation charges of the topological excitations of the two FFLO/PDW 
phases (rii and U are integers). 



Note that the expressions in Table |TT] are the most general expressions allowed by symmetry. For the GLW theory of 
Eq. [3 the elastic coefficients simplify. For example, for the phase ^ v - av , the elastic coefficients in the GLW limit are: 
Ps = 2|* | V, A = -(^) 2 |* | 2 K2 /2, n = (^) a |*o| a «i/2, 7 = (f ) 2 |*o| 2 («i + «a)/2, and e = ^|* | 2 K2 where * 
is the magnitude of the normalized order parameter. In general, higher order terms in the GLW theory will lead to 
new coefficients, for example the term 

w[(A q1 A Q2 )*(VAqi • VA Q2 ) + (A q1 Aq 3 )*(VA q1 ■ VA Q3 ) + (A Q iA Q3 )*(VA Q1 • VA Q3 ) + c.c] (6) 

will increase p s and decrease p and 7 (for positive w), leading to different energy scales for superfluid phase fluctuations 
and for phonons. Since terms such as Eq. [S] are not necessarily small near any melting transition in 2D, we have 
included all terms allowed by symmetry. We note that the related elastic theories for LO phases have been worked 
out microscopically in Refs. [lH and[2ll 

The elastic Hamiltonians in Table [TT| imply power law spatial correlations in 2D for the order parameters. While 
it is possible to carry out a complete analysis of the melting transition for ^a, this is not the case for fy v _ av . We 
therefore consider a simpler and more physically transparent approach that allows both "Fa and ty v - av to be treated 
on an equal footing. In particular, we ignore terms that give rise to spatial anisotropy in the correlation functions. 
This approach is akin to an early treatment of 2D melting done by Nelson^. This yields the correct phase melting 
diagram in 2D but does not provide accurate critical exponent a 23 i 24 . The simplified elastic Hamiltonian is 

H = f (W) 2 + ^n{Vu x f + (V % ) 2 ]. (7) 

The spatial dependence of the correlation functions is then given as 

(A Q (r)A_ Q *(0)) oc r -^+^ 

(Pq(r)P-Q(0)> oc r""- 
(p K (r)p- K (0)) « r" 3 ^ 
(/ z ,K(r)/ z ,- K (0)) cx r" 3 ™ 

(V 6e (r)V> 6 * e (0)) oc r" 9 ^ (8) 

where 77^ = T/{2irp s ) and rjd — T/(2np). 



IV. TOPOLOGICAL EXCITATIONS 



The low energy topological excitations are important in the melting of the FFLO/PDW phases. These are found 
by requiring single valuedness of the order parameter components Aq. In particular, defining Aq ; = Aoe 4 with 
6i = 9 + Q i ■ it, then along a contour surrounding a point defect in 2D, 

j> dO, = j> dO + <j> Q l ■ du = n^n (9) 

with integer rii. Implementing this condition for all components of the FFLO order parameter leads to the defect 
classification of Table [TTJ An important feature is the existence of defects that contain both fractional vortex charge 
and fractional dislocation charge. In the case of the stripe-like PDW and FFLO phases, the predicted 1/2 vortices 
play a central role in determining the phase diagra m 10 i 15 i 17 . Fractional vortices exist for ^ v ^ av but not for ^ a- 

The ground state supports: conventional vortices [n\ = n 2 = in Table HI)) : conventional dislocations 

/l(a,0) + Z 2 (f,^) (when m +n 2 +n 3 = in Table HI)) : and 1/3 vortices which combine a phase winding of 27r/3 and 
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a fractional dislocation (with a charge that is l/\/3 of the charge of the smallest conventional dislocation). We note 
that fractional vortices related to the 1/3 vortices found here have been discussed in the context of antiferromagnetic 
model on the XY lattice 2 ^, and the fully frustrated XY model on the dice lattice 2 ^. 

The ground state allows only conventional vortices and conventional dislocations. This latter result is somewhat 
surprising given the prevalence of fractional vortices in other FFLO/PDW ground states and is a direct consequence 
of the co-existence of the spatially uniform charge two s-wave superfluid with the FFLO/PDW order. 



V. DEFECT DRIVEN FFLO/PDW MELTING 



The elastic Hamiltonian of Eq. [7] implies that the interaction between the defects is given by: 



H top = 2TT^\K s n i n j + • 6,} ln (— ) (10) 



where m is the vortex charge of the defect, bi is the dislocation charge, a c is a hardcore cutoff, is the dis- 
tance between the defects i and j, K s = p s /T, and Kd = n/T. Core energies of the defects give rise to bare 
fugacities yi — exp{—C(K s n 2 + Kdbf/a 2 )} where C is a constant of order one. We consider the small fugac- 
ity limit. In describing the critical properties, only the lowest energy defects are required. For the ground state 
we include single vortices <f> dO = ±2ir (with fugacity y v ), minimal normal dislocations § du = a(±l,0), 
§ du = ±a(l/2, %/3/2), § du = ±a(— 1/2, %/3/2) (each with with fugacity yd), and one-third vortices § dO = ±2ir/3 
and § du = ±a^=(0, 1), ±a^=(^, ±a^=(^, |) (each with fugacity 2/1/3)- While for we include single 
vortices and minimal normal dislocations. 

The approach of Kosterlitz and Thouless as generalized to vector Coulomb gases^ r 24 ' 27 ' 28 is used to determine the 
phase diagram. The hard core cutoff a c of the defects is increased to a c — a c e dl . Under this infinitesimal coarse 
graining, pairs of defects separated by a distance a c either annihilate if they have opposite charges; or they combine to 
form a new defect described by the sum of the charges. This procedure leads to renormalization group (RG) equations 
for the scale dependent fugacities y^ and interaction parameters K s and Kd- For ^ v _ avi we find 



dK: 



dl 

dKd 1 
dl 
dy v 
dl 



47r 3 (y 2 + iy 2 /3 ) 

27r 3 (3y 2 + y 2 /3 ) 
(2 - TrK s )y v 



^ = (2 - nK d )y d + 2ny 2 d + 2^ 2 /3 



dyi/ 



'3 



dl 



= 2/l/3 {[2- 7 r(^ + ^)]+47ry d }. (11) 



To gain insight into Eq. [TTJ it is useful consider initially the first RG equation (for K s ). This RG equation arises due 
to screening of single vortices. These vortices can be screened by other single vortices and by the one-third vortices. 
The interaction between single vortices is given by K s , while that between a single vortex and a one-third vortex is 
given by K s /3. This leads to the first RG equation (for K s ): the 1/3 in the new term y 2 y 3 /3 comes from a factor 

of (K v /3) 2 (due to the interaction between a vortex and a one-third vortex) and a factor of three stemming from 
the three different possible ways to screen a single vortex with one-third vortices. The second RG equation follows 
from a similar consideration for dislocations (these can be screened by other dislocations and by one-third vortices). 
The equations determining the fugacities follow from the usual consideration s 27 ' 29 together with the possibility of 
combining defects to create a new defect^ -24 ' 28 . For example, dislocations can be created by pairing either two other 
dislocations or by pairing two one-third vortices. These two processes lead to the terms 27ry^ + 27ry 2 ^ 3 in the fourth 
RG equation. 

For comparison, we write the RG equations for ^a (these are simply the RG equations of superfluidity^ 7 - and 2D 
melting in the isotropic limit 2 ^), 

6^ 3 y 2 



dl 
dl 
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FIG. 3. Phase diagrams for and ty v -av- Phase boundaries are given by solid lines. The shaded triangles are discussed in the 
text. The new phases that arise from melting the FFLO/PDW phases are the density wave, superfluid, orbital antiferromagnet, 
and homogeneous charge six superfluid (SE^e) phases. The density wave phase has no superfluid order and breaks translational 
symmetry through the formation of a hexagonal lattice. The superfluid phase is a homogenous conventional paired superfluid. 
The orbital antiferromagnet phase has no superfluid order and hexagonal orbital current antiferromagnetic order. The *&6e 
phase is a spatially homogenous superfluid phase of bound states of six fermions 



^= (2-nK s )y v 

^ =(2 - nK d )y d + 2iry 2 d . (12) 

Despite starting with the same elastic Hamiltonian, the RG equations and the resultant phase diagrams for the two 
different ground states (shown in Fig. [3]) differ substantially. This is because the state ^ v - av has exotic one-third 
vortices, while the state "Fa does not. 

For the ground state \t a> the theory is that of uncoupled dislocations and vortices. Vortices proliferate for np s /T < 2 
and dislocations proliferate for irfi/T < 2. As a consequence, the melting of the FFLO/PDW phase occurs generally 
through two separate phase transitions. As temperature is increased, the first transition is either to a conventional 
superfluid or to a density wave state and the second is to the disordered phase. An examination of the correlation 
functions in Eq. [3] reveal that in the FFLO/PDW phase either the conventional superfluid order or the density 
wave order (and not the FFLO/PDW order) have the longest range correlations, masking the original FFLO/PDW 
order. Indeed, the melting of the FFLO/PDW to a density wave state and then to a disordered state represents the 
same sequence of transitions expected for a conventional supersoli d 30 i 31 . These arguments indicate that strongly 
resembles a conventional supersolid state. 

Let us come back to the state ty v - av . Its phase diagram can be understood qualitatively by considering the terms 
of Eq. [TT]that ar e linear in the fugacities. In this limit, vortices proliferate for irp s /T < 2 (superfluid order is lost), 
dislocations proliferate for np/T < 2 (density wave order is lost), and 1 /3 vortices proliferate when n(p s /9+fj,/3)/T < 2 
(all order is lost). Consequently, in addition to the disordered phase and the fully ordered FFLO/PDW phase, there 
exist two new phases. The first is an orbital antiferromagnet phase in which there is no superfluid order, however, 
the density wave order and the orbital order still exhibit power law correlations. The second phase is a charge 
six superfluid (^ee) that is spatially homogeneous. In this phase there is no density wave or orbital order. The 
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qualitative phase diagram for is shown in Fig. [3] (an exact phase diagram requires going beyond the small 

fugacity limit). The solid lines are the anticipated phases boundaries. These follow from an argument given in Ref.l32l 
for the phase diagram of 2D spinor condensates. In the shaded triangular regions of Fig. [3j the phases 4>6e and the 
orbital antiferromagnetic phase cannot be stable (a theory linear in fugacities would lead to the opposite conclusion). 
In particular, even though 7r(p s /9+/i/3)/T < 2 for the unrenormalized stiffnesses, this is not true for the renormalized 
stiffnesses. Consequently, 1/3 vortices will proliferate and all quasi-long range order will be lost. 

We note that, in principle, there can exist other phases that have not been explicitly considered here. For example, 
for rotationally invariant systems, the passage of the density wave phase into the disordered phase can have an 
intervening hexatic phased. Similarly, a recent analysis for the LO phase for rotationally invariant systems, leads to 
an intervening nematic phased. The consideration of such phases requires the inclusion of disclinations and related 
defects which have not been included in this work. 



VI. CONCLUSIONS 



We have presented an analysis of thermal melting on two FFLO /PDW ground states with hexagonal symmetry. For 
a FFLO/PDW vortex-antivortex lattice phase, we find that thermal melting can lead to either a charge six superfluid, 
an orbital antiferromagnetic, or directly to a disordered phase. While for a FFLO/PDW phase with only superfluid 
density oscillations, thermal melting necessarily proceeds in two transitions: the first to either a density wave phase or 
a conventional superfluid, and the second transition to a disordered phase. The latter FFLO/PDW phase is difficult 
to distinguish from a conventional supersolid phase. 

We acknowledge Manfred Sigrist for useful discussions. This work is supported by NSF grant DMR-0906655 and 
by Grants- in- Aid for Scientific Research (No. 19052003), MEXT of Japan. DFA and HT acknowledge the Hospitality 
of Pauli Center of the ETH-Zurich. 
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